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Gene-knockout experiments on single-cell organisms have estab- 
lished that expression of a substantial fraction of genes is not 
needed for optimal growth. This problem acquired a new dimen- 
sion with the recent discovery that environmental and genetic 
perturbations of the bacterium Escherichia coii are followed by 
the temporary activation of a large number of latent metabolic 
pathways, which suggests the hypothesis that temporarily acti- 
vated reactions impact growth and hence facilitate adaptation in 
the presence of perturbations. Here we test this hypothesis com- 
putationally and find, surprisingly, that the availability of latent 
pathways consistently offers no growth advantage, and tends in 
fact to inhibit growth after genetic perturbations. This is shown to 
be true even for latent pathways with a known function in alternate 
conditions, thus extending the significance of this adverse effect 
beyond apparently nonessential genes. These findings raise the 
possibility that latent pathway activation is in fact derivative of 
another, potentially suboptimal, adaptive response. 

complex networks | flux balance analysis | metabolic networks | gene dis- 
pensability I synthetic rescues 

Living cells are suqirisingly robust against mutations and, in par- 
ticular, against gene knockouts [1, 2, 3, 4, 5]. The origin of 
mutational robustness — whether it is a directly evolved trait or a 
byproduct of evolutionary history — remains debatable [6], In either 
case, metabolic network analysis shows that the nonessentiality of en- 
zymes and associated genes is largely due to the inactivity of the cor- 
responding metabolic reactions under laboratory conditions [7, 8, 9]. 
This leaves environmental robustness as the natural candidate to ex- 
plain gene nonessentiality. Yet, apart from chemical stress-based as- 
says [10], studies designed to test whether nonessential genes become 
essential under different conditions have failed to identify a pheno- 
type for more than a small fraction of additional genes [11]. A recent 
groundbreaking study has shown, however, that a large fraction of re- 
actions not active under standard laboratory conditions become tran- 
siently active after a genetic or environmental perturbation [12, 13]. 
Why? The prevailing interpretation has been that the transient activa- 
tion of such latent pathways facilitates adaptation to new conditions, 
thereby attributing function to genes that have been classified as dis- 
pensable for the lack of phenotype in steady-state experiments. This 
is naturally formulated as the hypothesis that latent pathways have 
a positive impact on postperturbation growth (cellular reproduction), 
which is a measure of competitive advantage with a strong empiri- 
cal basis [1, 2, 3, 10, 11, 13]. Even for genes with known functions 
under different conditions, this hypothesis is appealing as it suggests 
the possibility of an alternate phenotype that would not be detected 
in traditional high-throughput screens of knockout mutants [14]. 

Here we test this hypothesis using the most complete in silico 
reconstruction of the metabolic network of Escherichia coli K-12 
MG1655 [15, 16] and perturbations caused by single-gene knock- 
outs. The response of the metabolic network to knockout pertur- 
bations is modeled using both model-independent analysis and the 
two most accepted phenomenological approaches, minimization of 
metabolic adjustment (MOMA) [17] and regulatory on/off minimiza- 
tion (ROOM) [18] (Materials and Methods). Starting from a growth- 
maximizing state determined by flux balance analysis (FBA) [19], 



we compare the early postperturbation growth rate (Materials and 
Methods) of the original organism with that of a modified organism 
in which the latent reactions have been disabled. We consider glu- 
cose minimal medium and gene knockouts that necessarily change 
the original metabolic flux distribution but that nonetheless are com- 
patible with nonzero growth according to FBA. There are 52 enzyme- 
coding genes associated with 97 metabolic reactions in the recon- 
structed network that satisfy this condition. We systematically pre- 
dict the impact of latent pathway activation on growth rate following 
perturbations caused by the knockouts of each of these genes. 



Results 

Phenomenological analysis. Figure 1 illustrates the essence of our 
approach. In an optimal growth state, as observed for E. coli af- 
ter adaptive evolution in fixed environmental conditions [20], many 
metabolic reactions are inactive [8, 12]. Shortly after a perturbation, 
however, both the original and modified strain operate in a subopti- 
mal growth state [17, 18], which we model using MOMA or ROOM 
(Fig. lA). In silico [8] and laboratory [12] experiments show that 
this change is accompanied by a burst of reaction activity (Fig. IS), 
reflecting regulatory changes that locally reroute fluxes in the short- 
term metabolic response to the perturbation [21, 22]. If the per- 
turbation is nonlethal, the perturbed organisms will undergo adap- 
tive evolution — adopting beneficial mutations over longer timescales 
[23, 24] to achieve a new optimal growth state, which can be pre- 
dicted by FBA [25, 18]. 

For the perturbations considered in this study, the average and 
standard deviation of the number of transiently active reactions is 
291 ± 83 and 120 ± 59 for MOMA and ROOM, respectively. This 
difference is expected since ROOM, by design, favors a small number 
of significant flux changes, which reflects the fact that ROOM may 
model a later stage in the adaptive response pictured in Fig. IB than 
MOMA [18]. These numbers should be compared with the number 
of active reactions in the corresponding growth-maximizing states, 
which is 385 in the wild type and remains 385 on average in the 
knockout mutants, with an average of ~ 98% overlap between the 
two sets for the simplex solutions we consider; these numbers are 
representative for other choices of optima [26, 27] within our mod- 
els (Supporting Information, Sensitivity to Alternate Optima section). 
We emphasize that since the modified organism lacks only transiently 
active (latent) metabolic reactions, the optimal steady states are iden- 
tical to those of the unmodified strain both before and after the per- 
turbation. Our question is then whether the early postperturbation 
growth rate (before adaptive evolution) will be lower (red), remain 
equal (green), or become higher (blue) when these latent pathways 
are not present (Fig. lA). 

Our principal result is that the strains lacking latent pathways sys- 
tematically show equal or better adaptation as determined by growth 
within our in silico model, regardless of the approach used to model 
the organisms' response to perturbations. We assume that the organ- 
isms are in an optimal growth state both before and long after the 
perturbation, which accounts for cases that have received much at- 
tention in the literature [8, 12, 13], but we note that our conclusions 
remain equally valid when this assumption is relaxed (Supporting In- 
formation, Ejfects of Nonoptimal Reference States section). Table 
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Fig. 1 : Hypothetical growth impact of latent pathways under a knockout pertur- 
bation. (A), (6) The solid lines indicate the drop in growth rate {A) and the burst 
in latent reaction activation (6) that follow a gene knockout. The dashed lines 
indicate the possible behavior for a modified strain in which the original latent 
pathways have been removed: the postperturbation growth rate may decrease 
(red), remain the same (green) or increase (blue), and a smaller number of new 
latent pathways may be created. If the postperturbation growth rate is nonzero, 
the mutant is viable, and after a period of adaptive evolution it will converge to a 
new optimal growth state. This postperturbation optimal growth state is identical 
for both the modified and unmodified strain, and is characterized by a reduced 
number of active reactions relative to the suboptimal states. 



1 summarizes the results for all 52 single-gene knockouts consid- 
ered in our study. Relative to the unmodified strain in the subop- 
timal state following a gene knockout, the strain lacking the latent 
pathways exhibits equal or improved growth in 100% of the cases 
according to MOMA (Fig. 2 A) and in 98% of the cases according 
to ROOM (Fig. 2B). Across all knockout perturbations, this corre- 
sponds to an average change of +8.5% and +1.2% of the optimal 
wild-type growth rate, respectively. With either approach, a large 
fraction of mutants (50% for MOMA, 77% for ROOM) show a neg- 
ligible difference in growth rate (within ±1% of the wild-type growth 
rate) when the latent pathways are disabled. If only cases exhibiting 
significant changes are considered, the removal of latent pathways 



consistently increases the early postperturbation growth rate for all 
mutants, by an average of +16.9% for MOMA (Fig. 2A) and +4.6% 
for ROOM (Fig. 2B). Thus, for almost all knockouts, the strain lack- 
ing latent pathways is predicted to suffer no competitive disadvantage 
compared to the latent pathway-enabled strain. On the contrary, we 
predict that it more often shows improved growth in the suboptimal 
regime shortly after the perturbation. 

The set of transiently active (latent) reactions depends on the per- 
turbation. Even though we predict that in general, the removal of one 
of these 52 sets increases growth under the corresponding knockout, 
the same removal may in principle have an adverse effect under a 
different knockout. To address this possibility, we first note that the 
sets of simultaneously nonessential latent reactions remain sizeable: 
an average of 258 ± 79 for MOMA and 109 ± 53 for ROOM. For 
a given knockout perturbation, this set is defined as the subset of 
the original latent reactions that are inactive in the optimal growth 
state we consider for each of the other 51 knockout mutants. These 
reactions are therefore dispensable for optimal growth, both in the 
wild type and all 52 single gene knockout strains we consider, but 
are nonetheless transiently activated in response to the given pertur- 
bation. We have tested the impact of disabling these reduced sets 
of latent reactions under the corresponding knockouts (Materials and 
Methods). As shown in Fig. 2 C and D, the presence of these simulta- 
neously nonessential latent reactions has the same trend of inhibiting 
growth adaptation as found for the full sets of transiently activated 
reactions. 

The possibility that latent pathway activation enhances cells' via- 
bility following a perturbation is a compelling hypothesis, as it would 
reveal functions for genes that have thus far eluded high-throughput 
phenotype screens. We note, however, that our analysis also predicts 
the transient activation of pathways that do, in fact, have known phe- 
notypes under different conditions. For example, activation of the 
glyoxylate shunt is known to mitigate growth defects of E. coli on 
glucose following phosphofructokinase mutations [28]. Since we fo- 
cus on single knockouts, the genes affected by such mutations, pfliA 
and pfkB, are not among the perturbations we consider. Nonethe- 
less, out of the 52 unrelated knockout perturbations in our study, 
our models show the transient activation of the glyoxylate shunt in 
response to 25 of them according to MOMA and 7 according to 
ROOM. The same phenomenon can be observed for reactions that 
are essential under different environmental conditions but inactive in 
the aerobic glucose medium employed in our simulations. Pyruvate 
formate lyase is required for anaerobic growth in xylose medium ac- 
cording to experiments [29] and our models, but is transiently active 
for 2 (MOMA) and 18 (ROOM) of the 52 genetic perturbations in 
this study. This interesting effect — the nonspecific use of pathways 
under an array of perturbations quite different from the conditions 



Table 1: Summary of the predicted impact of latent pathways in £. co// K-12 MG1655 under single-gene knockout perturbations. 

MOMA ROOM Random 



Latent reactions for individual perturbations: 








All knockout perturbations 


-f8.5 (12.5)% 


-^1.2 (2.8)% 


(10.8)% 


Significant differences* 


-^16.9 (13.2)% 


-^4.6 (4.5)% 


+70.0 (10.8)% 


Number of reactions removed 


291 (83) 


120 (59) 


1,019 (5) 


Simultaneously nonessential latent reactions: 








All knockout perturbations 


+7 A (10.9)% 


-^1.2 (2.8)% 




Significant differences^ 


-^14.8 (11.4)% 


+5.2 (4.6)% 




Number of reactions removed 


258 (79) 


109 (53) 





* By more than 1 % of the wild-type growth rate: 50% (MOMA), 23% (ROOM), and 1 00% (random) of the perturbations. 
tBy more than 1% of the wild-type growth rate: 50% (MOMA), and 19% (ROOM) of the perturbations. 



Each column corresponds to the average (and standard deviation) of the difference in growth rate between the latent pathway-disabled and wild-type 
organisms for 52 different single-gene knockouts. The differences are expressed as percentages of the optimal wild-type growth rate. For all cases, 
the average postperturbation growth rate is higher for the strain without latent pathways. 
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Fig. 2: Predicted adaptive impact of latent pathways in £ co// under single-gene l<noci^out perturbations. Eacii row indicates tlie difference in early postperturbation 
growtti rate between tiie strains witti and wittiout latent pattiways for the l^nockout perturbation indicated along the vertical axis, when: {A and S) the modified 
organism lacks all reactions that are transiently active under the corresponding knockout perturbation; (Cand D) the modified organism lacks the set of simultaneously 
nonessential latent reactions, which is the subset of such latent reactions that are not needed for growth in any of the other 51 single-gene knockout mutants. For 
both MOMA {A and Q and ROOM (B and D), only cases showing a significant change in growth rate (> 1% of the wild-type optimal growth rate) are shown. Arrows 
indicate the extent of the increase in growth rate when the latent pathways are removed. The shaded background indicates the theoretical maximum growth rate for 
the mutant strain predicted by FBA. All growth rates are normalized by the optimal wild-type growth rate. The statistics are summarized in the left and center columns 
of Table 1 for MOMA and ROOM, respectively 



under which they have an observed phenotype — indicates that the 
phenomenon of latent pathway activation extends beyond the set of 
apparently nonessential genes. 

Model-independent analysis. The analysis above shows that the 
availability of latent pathways inhibits growth in the short term af- 
ter a genetic perturbation. But how sensitive are these conclusions 
to the models we used to simulate the response of the network? To 
provide model-independent evidence, we have determined how the 
volume of the space of feasible metabolic states (Materials and Meth- 
ods) depends on the growth rate. As shown in Fig. 3 A and B (green 
lines) for the cyoA and Ipd knockout mutants, the volume system- 
atically decreases as a function of growth rate for the single-gene 
knockout mutants considered in our study, indicating that the number 
of metabolic states available to the unevolved mutant is much larger 
at lower growth rates. When the latent reactions are disabled, how- 
ever, the relative volume, and hence the relative number of available 
metabolic states, increases for large growth rates (Fig. 3 A and B\ blue 
lines). Therefore, the principal effect of removing latent pathways ap- 
pears to be an increase in the relative frequency of high-growth states 
due to the preferential elimination of low-growth states. It should 
be noted, however, that a large number of high-growth states are 
also disabled in this process due to the "entanglement" between la- 
tent pathways and biomass-producing pathways that exist under the 
metabolic steady-state conditions of our models (Supporting Infor- 
mation, Elementary Mode Analysis section). 

To appreciate the constraints imposed by this structure, imagine 
that the organism responds to perturbations by moving to a random 
metabolic state in the feasible space of fluxes. We simulated this hy- 
pothetical response using an implementation of the hit-and-run sam- 



pling algorithm {Materials and Methods). As shown in Fig. 3C, the 
postperturbation growth rate is nearly zero for all mutants with latent 
pathways and close to the theoretical maximum for all mutants with- 
out them. This random response is arguably a lower bound for the 
actual response of organisms that have evolved to cope with pertur- 
bations, but the conclusion is clear: unless we assume that organisms 
have evolved to respond to perturbations in a highly specific manner, 
which appears to be inconsistent with experiments [30], the availabil- 
ity of latent pathways does not facilitate growth, and this prediction 
is largely independent of the network response to perturbation. This 
holds true in particular for MOMA and ROOM, which incorporate 
(in different ways) the main flux rerouting features observed in the 
activation of latent pathways in E. coli [12]. 

Further mechanistic insight comes from the recently identified 
synthetic rescue interactions [31], in which the knockout of a gene 
inhibits growth but, counterintuitively, the targeted concurrent knock- 
out of additional genes recovers the ability of the organisms to grow. 
The reactions catalyzed by such rescue genes are predicted to be ac- 
tive in typical suboptimal states and inactive in growth-maximizing 
states of the knockout mutant [8, 31]. Now, given the observation 
above that the set of active reactions predicted by FBA is only slightly 
modified by a gene knockout, it follows that most rescue genes are 
in fact associated with latent pathways. This in turn explains why 
the latent pathway-disabled strains show improved growth. Note that 
this argument cannot be anticipated from intuition, because an enor- 
mous number of low-growth states (up to several orders of magnitude 
larger than for near-optimal growth) may exist even when latent path- 
ways are disabled (Fig. 3 A and B). Furthermore, even in the extreme 
case when one disables all reactions that are inactive in the optimal 
state of the knockout mutant, metabolic states with a very low growth 
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Fig. 3: Properties of the space of feasible metabolic states. {A and 6) Volume of the metabolic solution space as a function of the growth rate for the cyoA- (A) and the 
Ipd- (6) knockout mutant. All values are normalized so that the area under each curve is 1 . Thus, each curve can be interpreted as a probability distribution of growth 
rate within the corresponding metabolic space, either with or without latent pathways. For increasing growth rate, the normalized volume decreases systematically 
when the latent pathways are available (green lines), but the curve becomes less steep and the volume may even peak at nonzero growth when latent pathways 
are disabled (blue lines). These results were obtained for the MOMA-predicted response of the network, and similar behavior is observed for the other single-gene 
knockout mutants in our model. Disabling latent pathways therefore increases the relative frequency of high-growth metabolic states at the expense of low-growth 
states. For numerical feasibility, this calculation was implemented using the central metabolism of £ coli {Materials and Methods). (Q Postperturbation growth rate 
if the network were to respond by moving to a random position in the space of feasible states. The dots and error bars correspond to the average and standard 
deviation for each of the knockout mutants indicated at the bottom, when latent reactions are available (green) and when latent pathways associated with the particular 
perturbation and random response are disabled (blue). Interpreting a random response as a lower bound for the likely response of the organisms, the systematically 
higher growth of the latent pathway-disabled strains corroborates the conclusion that latent pathways do not facilitate adaptation. These statistics are summarized in 
the right column of Table 1 . All growth rates are normalized by the FBA-predicted rate of the knockout mutants. 



rate (« 10% or less of the wild-type optimum) exist in 47 out of the 
52 mutants (Supporting Information, Elementary Mode Analysis sec- 
tion). This threshold is significant since all but 5 unmodified knock- 
out mutants exceed this growth rate according to MOMA, and all but 
1 according to ROOM. Thus, although our model-independent anal- 
ysis suggests that latent pathway activation inhibits growth under a 
general choice of metabolic state after a perturbation, this is not di- 
rectly imposed by the geometry of the solution space. Rather, the 
predicted growth benefit associated with latent pathway removal and 
synthetic rescues is partly a reflection of the cells' adaptive response. 

An extreme example of this rescue effect is provided by the cyoA- 
deficient mutant, which is predicted by MOMA to drop to < 10% of 
the optimal wild- type growth rate following the perturbation, but re- 
covers to ~ 60% if the latent pathways are also disabled. In addition, 
cases in which the single-gene knockout mutant operates near the 
theoretical optimum but growth nonetheless improves upon the re- 



moval of latent pathways, such as for the/o/D mutant, can be related 
to weaker foms of synthetic rescues [31]. 

This surprising relation to synthetic rescues is particularly in- 
teresting when we note that latent reactions define several pathways 
whose participation in E. coifs metabolism has been controversial or 
elusive. The Entner-Doudoroff (ED) pathway, an alternative to gly- 
colysis for glucose catabolism, is inactive in wild-type E. coli accord- 
ing to in vivo experiments in glucose [32] but becomes transiently 
active in mutants lacking phosphoglucose isomerase [12]. This ac- 
tivation may serve to reduce NADPH accumulation accompanying 
increased flux through the pentose phosphate pathway [33]. Both 
MOMA and ROOM predict a small, transient flux through the ED 
pathway in response to the knockout of pgi, the gene coding phos- 
phoglucose isomerase. In triphosphate isomerase-deficient strains, 
our model predicts the activation of the normally inactive methyl- 
glyoxal bypass [34]. Experimentally, this pathway is observed to 
channel excess dihydroxyacetone phosphate (DHAP) into pyruvate 
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following glycolytic flux splitting into glyccraldehyde 3-phosphate 
and DHAP after the knockout of the associated gene, tpi [12]. These 
findings emphasize the importance of probing multiple gene knock- 
outs or perturbations — previously suggested in the context of syn- 
thetic lethality [35], synthetic rescues [31, 36], multi-target drug dis- 
covery [37, 36], and neutral mutations [38] — as a means to determine 
the puzzling role of transients. 



Discussion 

Latent pathways and their associated genes are, by definition, dis- 
pensable for optimal growth under the given conditions both before 
and after a perturbation. Several explanations have been offered to 
justify the persistence of apparently nonessential metabolic genes 
within the genome. Environmental robustness is the most natural 
hypothesis, as it acknowledges the unpredictable, time-varying con- 
ditions that confront single-cell organisms in natural environments. 
Yet, in silico [7, 39] and experimental [11] studies of model organ- 
isms under various environmental conditions likely to occur in nature 
suggest that condition-dependent robustness is inadequate to fully ex- 
plain the dispensability of metabolic genes. Part of the remaining re- 
dundancy may be attributed to the varying efficiencies of redundant 
pathways under different environments [40], or selective pressure for 
increased metabolic capacity across all conditions [39]. An alterna- 
tive to environmental robustness is genetic robustness, in which re- 
dundant pathways buffer against null mutations [41, 42]. But regard- 
less of which explanation might apply in the context of latent path- 
ways, the question that follows from our analysis is not why these dis- 
pensable pathways are present in the genome, but rather what causes 
their transient activation in response to genetic perturbations. Our 
results indicate that this activation confers no advantage in fitness as 
measured by growth, and more often hinders growth in the short term 
following a perturbation. 

The dispensability of latent pathway activation predicted in this 
study can be interpreted in three different ways, which are not neces- 
sarily mutually exclusive. First, it is possible that temporary reaction 
activation does not provide necessary intermediate states but is in- 
stead a byproduct of the network's suboptimal response to perturba- 
tions. The fact that E. coli undergoes a period of suboptimal growth 
following genetic perturbations is well-established in both in silico 
and in vivo experiments [17, 18]. This interpretation is thus sup- 
ported by the recent observation that typical suboptimal metabolic 
states have a much larger number of active reactions (more than 2.5 
times larger for the conditions considered here) than typical states 
that optimize growth [8]. 

A second interpretation is that cells activate a large number of 
reactions as a global nonspecific response to perturbations that nev- 
ertheless creates a library of possible metabolic states that can be 
subsequently fine-tuned by adaptive evolution. This scenario can- 
not be rejected by existing experimental results [30] and is appealing 
as it allows for an indirect, long-term benefit of temporary reaction 
activation even if, as predicted here, this activation inhibits growth 
in the short term. Although whole-cell regulation remains largely 
unexplained, there are known mechanisms that could subsequently 
lead to optimal growth [43, 44]. This interpretation suggests that 
fine-tuning of whole-cell reaction activity is best achieved by down- 
regulating specific over-expressed reactions rather than by the coor- 
dinated up-regulation of entire under-expressed pathways, a principle 
that is widely appreciated in metabolic engineering but that is yet to 
be demonstrated for natural systems. The recent observation that E. 
coli's response to perturbations is more stable from the view point of 
metabolite concentration and protein or mRNA abundance [45] than 
reaction flux [12] may be relevant for the validation of this interpre- 
tation. The objection to the second interpretation is that it cannot 
explain (not even in the long term) the availability of pathways that 
are simultaneously dispensable for all known pertiu-bations that have 
shaped the evolution of the organisms. This "paradox of latency" re- 



inforces the need to study the mechanisms that govern latent pathway 
activation. 

In experimental studies of gene function in microbial organisms, 
growth is the most often used indicator of fitness, owing to both its 
accessibility for measurement and its intrinsic importance in deter- 
mining viability [1, 2, 3, 46]. Although the results of this study sug- 
gest that the availability of latent pathways does not promote faster 
growth following a perturbation, the conclusion that this negatively 
impacts adaptation is reached by studying these pathways in the con- 
text of metabolism alone. Thus, a third interpretation is that, in gen- 
erd, the cellular objective invoked in the adaptive response is not 
growth. Instead, the activation of otherwise latent metabolic path- 
ways may accompany other cellular processes, either inside or out- 
side metabolism, that are initiated to ensure survival. For example, 
yeast adopts changes in cell shape and internal pressure in response to 
osmotic shock, a process that recruits the metabolism to accumulate 
specific metabolites [47]. Moreover, enzymes and metabolic cofac- 
tors involved in transiently upregulated pathways may have regula- 
tory or signaling roles in addition to their metaboUc function [48]. 
This third possibility allows for an advantageous, external function 
for pathways whose activation appears disadvantageous when only 
the metaboUsm is considered. In this way, latent pathway activation is 
incorporated into a larger, more sophisticated adaptive response. This 
explanation, however, seems inconsistent with our observation of the 
apparently highly nonspecific activation of pathways in response to 
different perturbations. 

Regardless of which interpretation proves to be correct, this study 
leads to a clear, experimentally testable prediction: that latent path- 
way activation does not enhance, and in fact often inhibits, early 
postperturbation growth. While we expect this behavior to be ob- 
served experimentally under a wide range of conditions, regardless 
of the specific suboptimal response of the cell, deviations from this 
behavior would also be highly informative since they would uncover 
growth phenotypes not detected in previous steady-state experiments. 
Varying environments can also generate strong [45] and sometimes 
counterintuitive responses, such as the possibility of accelerating the 
evolution of unevolved strains [49]. Therefore, it may well be the 
case that in time-varying environments, which are beyond current 
modeling capability, latent pathway activation will exhibit a different 
fitness effect. 

Materials and Methods 

Network and perturbations. We used the /AF1 260 E. coli model [1 5], which Is 
the most up-to-date reconstruction of the metabolic network of E. CO//MG1655. 
The network consists of 2,082 unique reactions catalyzed by 1 ,260 genes and 
Involves 1,369 metabolites, as well as 299 exchange fluxes and the biomass 
flux. We focused on the 52 single-gene knockout strains that are compatible with 
growth but for which the original growth-maximizing flux state becomes Infeasl- 
bie after the gene knockout. The volume calculation, which remains challenging 
for the full network, was performed using a reduced network consisting of 62 
reactions, 101 genes, 49 metabolites, and 14 exchange fluxes representing the 
central metabolism of E. coli [50]. 

Medium and constraints. The simulated medium consisted of limited amount 
of glucose (8 mmoi/g DW-h) and oxygen (18.5 mmoi/g DW-h), and unlimited 
amount of carbon dioxide. Iron (II), protons, water, potassium, sodium, ammonia, 
phosphate, and sulfate. Irreversible reactions are constrained to have nonneg- 
atlve fluxes. The flux through the ATP maintenance reaction was set to 8.39 
mmol/g DW-h. A total of 1,432 reactions in the /AF1260 model, including 61 
reactions in the subnetwork of the central metabolism, can be active under these 
medium conditions. An analysis of the effect of regulatory limitations [51 ] that may 
constrain this reaction activity in vivo is considered In the Supporting Information 
{Regulatory Constraints section). There It Is shown that our results remain valid 
under the additional constraints Imposed by these limitations. 

Feasible metabolic states. The state of the metabolic network is represented 
by a vector of ail reaction fluxes v = (i/j). Since the time scale at which the 
network responds to perturbations is much shorter than the characteristic time 
for adaptive evolution, we focused on steady-state flux distributions both before 
and after perturbation. Steady-state fluxes are solutions of the mass-balance 
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equation S ■ v = 0, where S = (Sy) is the matrix of stoichiometric co- 
efficients, under the constraint imposed by the medium, reaction irreversibiiity, 
ATP maintenance requirements, and possibiy gene l<nocl<outs. A l<nocl<out of the 
enzyme-coding genes associated with reaction j is implemented through the ad- 
ditional constraint Vj = 0. The solution space is the set of all such steady-state 
flux vectors and it has the form of a convex polytope. We refer to the individual 
solutions in this space as feasible metabolic states. 

Objective functions. FBA [19] identifies a growth-maximizing state within the 
space of feasible metabolic states by maximizing the flux ut, through a biomass re- 
action that drains biomass precursors. With respect to this original state, MOMA 
[17] and ROOM [18] identify feasible states that minimize the distance in the 
space of fluxes and the number of significant flux changes, respectively (Sup- 
porting Information, sections 2-3). Our implementations of FBA, MOMA, ROOM, 
and the hit-and-run algorithm used a commercial optimization package (ILOG 
CPLEX, version 11.0, www.ilog.com). For all FBA results, we have used the 
growth-maximizing states provided by the simplex algorithm. More information 
about the computational methods is provided in Supporting Information, where 
it is also shown that our results do not depend sensitively on the assumption 
of optimal growth in the reference states either before or after the perturbation 
{Effects ofNonoptimal Reference States), nor on particular choices for the growth- 
maximizing states used throughout the paper {Sensitivity to Alternate Optima). 

Biomass flux and growth rates. The in silico model predicts the biomass 
flux, but for exponential growth the result can be expressed in terms of a normal- 
ized growth rate ii. This follows from the observation that biomass production is 
governed by ^ = k, where k is the growth rate, N measures the popu- 
lation size, and ^ is proportional to the biomass flux Vb. Therefore, when 
normalized with respect to the wild-type or theoretical maximum, the normalized 
biomass flux Vb/vb,o equals the corresponding normalized growth rate k/kq 
used throughout the paper. 

Identification of iatent pathways. We define the latent reactions (pathways) 
associated with the l^nockout of gene A to be the set La of all reactions predicted 
to be transiently active in the unevolved mutant shortly after the perturbation (ac- 
cording to MOMA, ROOM, or the hit-and-run algorithm), but inactive in both the 
optimal wild-type and mutant strains (according to FBA). This set is therefore 
nonessential for optimal growth under the knockout of A, although in principle 
some of these reactions may be necessary for growth under the knockout of a 
different gene, B. To account for this, we tested the impact of a set L'^ of latent 
reactions that are simultaneously nonessential in the optimal growth states of the 
other 51 single-gene knockout mutants that we consider. 



Voiume caicuiation. The exact calculation of the volume of the (high- 
dimensional) solution space is computationally intractable and, because this set 
is very skewed [52], even approximate calculations are computationally demand- 
ing. To determine the volume of the solution space as a function of growth rate, 
we used an approximate inference algorithm based on Belief Propagation [53]. 
In this approach, the convex polytope representing the constrained flux space is 
tiled with hypercubes of size e. We then use a message-passing algorithm to 
approximate the probability distribution P{i') over the discretized space. Using 
this, we define the associated entropy S = —^2 P{^) logio ^i'^)- w^iich 
counts the log of the number of e-hypercubes that overlap the space of feasible 
states. The total volume covered by these cubes is then V = lO'^g", where n 
is the dimension of the space. If m is the number of linearly independent mass 
balance constraints and / is the number of available fluxes after a given knockout, 
then the dimension of the unmodified metabolic space \sn = f — m while the 
solution space for the modified organism has dimension n = f — m — l, where 
/ is the number of independent latent reactions removed. We used a granularity 
of e = ^ for all calculations presented in the paper. 

Hit-and-run aigoritlim. To randomly sample the metabolic solution space, we 
implemented a hit-and-run algorithm [54] with artificial centering [55], which is 
a quickly converging sampler for high-dimensional convex sets. The algorithm 
is based on selecting a randomly oriented line I passing through the current 
sample point, and then selecting the point for the next iteration from a uniform 
distribution along I. For efficient sampling, we employed artificial centering [55], 
where the orientation of the line I is obtained from the direction defined by the 
current sample point and the center (mean) of a subset C of already-sampled 
points. The subset C was initially created by taking 10,000 warm-up points on 
the boundary of the space, and was updated recursively by replacing a randomly 
selected point of C with the currently sampled point. In all calculations presented 
in the paper, we generated a set of 5 x 1 0* points to sample the solution space. 
All calculations were performed with the COBRA Toolbox [56]. 
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Supporting Information 

A metabolic network of m metabolites and n reactions is represented 
by a m X n matrix S = {Sij), where Sij is the stoichiometric coeffi- 
cient of metabolite i in reaction j. The state of the metabolic network 
is represented by a vector of reaction fluxes i/ = (vj), where i/j is 
the flux through reaction j. In our analysis, the steady-state solutions 
of the system are determined by mass balance constraints, 

Si^ = 0, [SI] 

and additional constraints that limit the range of the individual fluxes, 

[S2] 

where the inequalities in this notation are assumed to apply to each 

component individually. The bounds on individual fluxes are deter- 
mined by substrate availability in the given medium, the ATP mainte- 
nance requirement, and thermodynamic constraints that limit the re- 
versibility of the corresponding reaction. The knockout of the genes 
associated with enzymes catalyzing reaction j corresponds to the ad- 
ditional constraint Vj = 0. The exchange fluxes and the biomass flux 
are excluded from our implementation of the MOMA and ROOM 
objective functions. 

Flux Balance Analysis (FBA). FBA [1] is used to identify flux vec- 
tors 1/ satisfying (SI) and (S2) that maximize biomass production, 
which is represented by an additional reaction b that drains biomass 
components. The problem is implemented as a linear program: 

max Ub [S3] 
s.t. 5 ■ 1/ = 

l^min ^ ^ l^max. 

In general, the FBA solution is not unique, and the results in this pa- 
per were obtained by selecting the optimal solution provided by the 
simplex algorithm. For a discussion of the effects of choosing alter- 
nate FBA solutions, see Sensitivity to Alternate Optima below. 

Minimization of lUletabolic Adjustment (IVIOIVIA). MOMA [2] selects 
a suboptimal growth state i/ by minimizing the Euclidean distance 
in flux space from the wild- type optimal growth state, w. This is 
implemented as a quadratic programming problem: 

min {v — w)^ ■ (v — w) [S4] 
s.t. Sv = 

Vj =0, j € A, 

where A is the set of indices corresponding to reactions deactivated 
by gene knockouts. 

Regulatory On/Off Minimization (ROOM). UnUke MOMA, which 
favors a potentially large number of small-magnitude flux changes, 
ROOM [3] chooses a suboptimal growth solution v with a minimal 
number of "significant" flux changes from the original state. This is 
usually implemented as a mixed-integer programming problem (in- 



teger ROOM): 

71 

min [S5] 
s.t. S-v = Q 

^min ^ ^ '^max 

i'k = 0, k e A 

forj = l,...,n 

toJ = Wj + 5\Wj \ + € 

w\ = Wj - 6\wj I - e 

yj e {0, 1}, 

where w and A are as for MOMA, and S and e express tolerances for 
relative and absolute change from the original state, respectively. 

In our numerical experiments we chose to use a linear program- 
ming variant of the method (linear ROOM) obtained by allowing the 
above binary constraints to be continuous in the interval < < 1 
and setting S = e = 0. Linear ROOM is biologically well-motivated, 
given that gene activity is best described by a continuous variable, 
and has the advantage of being computationally inexpensive for all 
mutants. Table SI shows a comparison of the growth impacts pre- 
dicted by integer ROOM and the linear variant upon latent pathway 
removal. For most mutants, the two methods show similar change 
in growth rate (increase, decrease, or an insignificant change) when 
the latent pathways are removed, even though the exact growth rate 
predictions may differ. 

For the integer variant, we followed Shlomi et al. [3] in choosing 
the values 5 = 0.03 and e = 0.001, which yielded reasonable run- 
ning times for most knockout strains. Even for larger values of these 
tolerances values there are a handful of cases for which no optimal 
solution is found in any reasonable amount of time. For these cases, 
in the comparison of Table SI, we take the best solution found after 
one hour of computation on a 3.4 GHz CPU. The integer ROOM solu- 
tions were further constrained to minimize the aggregate flux change 
from the wild-type optimal state, — Wj\. 

Sensitivity to Alternate Optima. In general, the optimal flux distribu- 
tion given by FBA is not unique, neither before nor after a given gene 
knockout [4, 5]. The implications of this nonuniqueness must be con- 
sidered for two reasons. First, the set of transiently active reactions is 
defined with respect to the exact optimal flux distribution before and 
after the knockout perturbation, as well as the suboptimal flux dis- 
tribution after the perturbation (see Effects of Nonoptimal Reference 
States below for an analysis of the effects of choosing suboptimal 
reference states in this definition of latent pathways). Second, the 
postperturbation flux distribution is itself dependent on a particular 
choice for the original optimal state, since both MOMA and ROOM 
operate by minimizing some distance to a reference flux state. 

To test the sensitivity of our results to a particular choice of FBA 
solutions, we repeated our simulations for numerous combinations 
of wild-type and optimal mutant flux distributions. Wc sampled the 
available FBA-predicted states by fixing the corresponding growth 
rate (either wild type or optimal mutant) and maximizing or min- 
imizing each of the reaction fluxes allowed to vary under this ad- 
ditional constraint. After choosing a particular pair of optimal flux 
distributions in this way, we determined the corresponding set of la- 
tent pathways according to MOMA and ROOM, and the associated 
suboptimal growth rates before and after the removal of these latent 
pathways. Figure SI shows the results of this analysis for the ppk- 
and fp;A-knockout perturbations. The distributions of growth rate 
predictions for the strains with (green) and without (blue) latent path- 
ways do not overlap. A less extensive sampling for each of the other 



8 



50 knockout mutants considered in our study shows similar behav- 
ior. Altogether, these results suggest that our predictions about the 
growth effect of latent pathway availabiUty are robust with respect to 
alternate optimal flux distributions. 

Regulatory Constraints. All results for the full E. coli iAF1260 
model presented in the main text were computed under the uptake 
and steady-state constraints listed in the Materials am! Methods sec- 
tion. These constraints do not take into account regulatory effects, 
which may limit the set of genes that can be transcribed under the 
nutrient conditions we consider. Regulatory effects are thus expected 
to constrain the set of steady-state metabolic states predicted by FBA, 
MOMA, and ROOM that can actually be realized in vivo [7]. To ad- 
dress the possible impact of regulatory constraints on our predictions, 
we have repeated all calculations in the paper for a modified version 
of the iAF1260 model, following Feist et al. [8] in disabling a set 
of 152 reactions beforehand. There is evidence that, due to regula- 
tory constraints, these reactions are inactive in the aerobic glucose 
conditions we simulate. Under these additional constraints, there are 
46 single-gene knockouts that change the original flux distribution 
but nonetheless allow nonzero growth in the resulting knockout mu- 
tants according to FBA, which is the same criterion we used to select 
knockout perturbations in the main text for the unmodified /AF1260 
model. 

Table S2 summarizes the predicted growth impact of disabling 
latent pathways in the modified /API 260 model. Regardless of the 
approach used to simulate the response (MOMA, ROOM, or hit-and- 
run sampling), all 46 knockout mutants in this model show nearly 
equal (within ±1% of the wild type) or improved growth in the short 
term following the perturbation when the latent pathways are dis- 
abled. Moreover, the average growth increases and the numbers of 
transiently activated reactions are comparable to those presented in 
Table 1 of the main text for the unmodified iAF1260 model. Our 
results are therefore not dependent on having the full complement 
of metabolic reactions in the network available for activation in the 
organisms' initial response to perturbations. Indeed, the predicted 
adverse growth effect of latent pathway activation is expected to hold 
even under additional constraints that may reflect other limitations of 
the metabolic response in vivo. 

Effects of Nonoptimal Reference States. The phenomenological 
models used in the main text predict a suboptimal metabolic state by 
minimizing the value of a distance metric with respect to the preper- 
turbation reference state, which is assumed to be growth-maximizing 
as predicted by FBA. The postperturbation reference state was also 
assumed to be growth-maximizing after adaptive evolution. These 
simplifications overlook the possibiUty that in many natural environ- 
ments the metabolic state may be nonoptimal even before and long 
after an external perturbation. This situation could arise, for example, 
under time- varying conditions, where environmental changes prevent 
the organisms from approaching optimal growth states. In this case, 
some of the pathways we have classified as latent may fail to qualify 
in vivo since they carry nonzero flux in one or both of the reference 
states. 

Geometrically, this scenario corresponds to the appropriate 
(nonoptimal) reference states lying in the interior of the feasible 
metabolic solution space, where many more reactions are active, as 
opposed to the boundary of the space, where reaction activity is lim- 
ited by irreversibiUty constraints [9]. This situation can be accom- 
modated with a variant of our modeling approach. We have system- 
atically tested the growth impact of reaction upregulation relative to 
nonoptimal reference states, which is the natural generalization of the 
latent reaction activation considered in our original analysis. Given 
two optimal states vi (wild type) and 1^2 (evolved mutant) and their 
associated biomass fluxes, v'^'^^ and fc'^j*' we now limit the biomass 
flux to some fraction of these optimal values. This is imposed within 



our in silico models by the additional constraints z/t i = Az/^^ and 

Vb,2 ~ Aj/j 2*' where < A < 1. The nonoptimal reference states 
are then defined by replacing the optimal states before and long after 
a knockout with the closest feasible metabolic states that satisfy these 
growth constraints. We performed this analysis for A — 0.4 and 0.7. 

With respect to the new choices of nonoptimal reference states, 
the short term response of the metabolic network to single-gene 
knockouts still exhibits a transient burst of reaction activity. Accord- 
ing to MOMA, the average and standard deviation of the number of 
fluxes with larger magnitude than in both reference states is 260 ± 83 
(A — 0.4) and 263 ± 82 (A = 0.7). These numbers are comparable 
to the 291 ± 83 reactions that are transiently activated by MOMA 
with respect to optimal reference states (main text, Table 1). Figure 
S2 shows the MOMA-predicted difference in growth rate when these 
transiently upregulated fluxes are constrained to not exceed the ref- 
erence states in magnitude. In all cases, downregulation of the tran- 
siently upregulated pathways is predicted to improve growth in the 
short term following a knockout perturbation, by an average of 6.0% 
(A = 0.4) and 10.4% (A = 0.7) of flie optimal wild-type growth 
rate. This analysis is therefore in agreement with the prediction pre- 
sented in the main text for optimal reference states, namely, that the 
transient activation (or otherwise upregulation) of latent pathways 
generally inhibits growth in the short term following a perturbation. 

Elementary Mode Analysis. The results of Fig. 3 (main text) sug- 
gest that the primary effect of latent pathway removal is to favor the 
availabihty of high-growth metabolic states by preferentially elim- 
inating low-growth states following a genetic perturbation. This is 
accomplished by eUminating reactions that are silent in the optimal 
states we consider before and after the knockout, thereby increasing 
the likelihood that the initial metabolic response will activate path- 
ways associated with higher growth. It should be noted, however, 
that a given reaction can in general be active in many metabolic states, 
spanning both low- and high-growth phenotypes. It is therefore likely 
that many high-growth states will by eliminated by disabling latent 
pathways as well. 

To systematically examine this connection between high-growth 
states and latent pathways, we have used Elementary Mode (EM) 
Analysis, which is an approach for analyzing metabolic networks in 
terms of interconnected sets of reactions [10]. An EM is defined as a 
unique set of active reactions in the network (represented by a vector 
1/ in flux space) that ;) satisfies the steady-state constraints S -v = 0, 
ii) obeys all reversibility constraints (negative entries in u must cor- 
respond to reversible reactions), and Hi) is minimal in the sense that 
no reaction may be removed from the set while still satisfying ( i) and 
(ii) [11]. Any steady-state flux distribution can be represented as a 
linear combination of EMs with noimegative coefficients. 

The number of EMs of grows combinatorially with the size of 
the metabolic network, making their calculation computationally in- 
feasible for the full E. coli /API 260 model. Therefore, we focused on 
the subnetwork comprising E. coli's central metabolism, which was 
also used for the volume calculation (main text, Model-Independent 
Analysis). Using the program METATOOL [12], we obtained the 
full set of 18,656 EMs available on glucose. We have classified each 
mode as "biomass-producing" or "nonbiomass-producing", based on 
whether it has a positive or zero entry corresponding to the (irre- 
versible) biomass flux, respectively. It follows that a linear combi- 
nation of EMs representing a general zero or low-growth metabolic 
state will be composed primarily of non biomass-producing EMs, 
with only a small aggregate contribution from the EMs that produce 
biomass. Figure S3 shows the effect of latent pathway removal on 
these two types of EMs for the cyoA and Ipd mutants. A significant 
fraction of the original 18,656 EMs are eliminated by the knockout 
perturbation (namely, those that involve a disabled reaction). Further 
EMs will be eliminated when the latent pathways are disabled. As 
expected from the solution space volume calculation and hit-and-run 
analysis in the main text {Model-Independent Analysis section), the 
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effect of latent pathway removal is to skew the distribution of EMs 
toward those that produce biomass (Fig. S3A). But, surprisingly, a 
large number of biomass-producing modes are sacrificed as well, and 
in fact comprise the majority of EMs disabled in this process (Fig. 
S3 A and B). 

Additional insight comes from analyzing the growth capabilities 
of the metabolism when one eliminates every reaction that is silent in 
the corresponding FBA-predicted optimal state of the knockout mu- 
tant. This corresponds to an average of 1,967 ± 6 reactions across 
all 52 knockout mutants, roughly 7 and 16 times larger than the num- 
ber of latent pathways removed for MOMA and ROOM, respectively 
(main text. Table 1). Figure S4 shows the range of growth rates that 
can be realized by each mutant in this scenario. For the majority of 
the mutants, metabolic states with very low growth (< 10% of the 
optimal wild-type growth rate) exist, even when this large set of op- 
timally silent reactions is disabled. This is particularly significant 
given that, with few exceptions (5 according to MOMA, 1 according 
to ROOM), the unadapted states of the 52 tested knockout mutants 
exceed this growth threshold in our models when the latent pathways 
are enabled. 
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Fig. S1: Impact of alternate optima. {A and 6) Sensitivity to alternate optima tor the pp/<-knockout perturbation according to MOMA (A) and ROOM (6) (based on 
12,760 different pairs of optimal states). (Cand D) Sensitivity to alternate optima for the tp//4-knockout perturbation according to MOMA (Q and ROOM (D) (based 
on 20,880 different pairs of optimal states). Each panel shows the probability distributions of the predicted early postperturbation growth rates for strains in which 
the latent pathways are available (green) and disabled (blue). All growth rates are normalized by the optimal growth rate of the wild type. Dashed lines indicate the 
growth rate predictions presented in the results of the main text. In all cases, the growth rates are tightly clustered around these values, and the two distributions do 
not overlap. 
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Fig. S2: Impact of nonoptimal reference states. {A and 6) Effect of latent pathway upregulation using reference metabolic states (wild type and "adapted" mutant) 
whose growth rates are limited to 70% {A) and 40% (6) of the respective optima. For each knocl^out, the panels show the change in the MOMA-predicted growth rate 
after the perturbation when we limit the magnitude of the flux through each latent reaction to the corresponding highest value between the two reference states. The 
dotted lines and gray bars indicate the maximum allowed growth rates for the wild type and knockout mutants, respectively. All values are normalized as a fraction of 
the maximum wild-type growth rate. We only display cases that show a significant difference in growth rate (> 1% of the wild-type optimum) when the latent pathways 
are downregulated. 
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Fig. S3: Effect of latent pathway removal on elementary modes of £. co//'s central metabolism. These simulations are based on the Ipd- and cyo/\-knockout mutants. 
{A) Number of elementary modes having a zero (gray) and positive (blue) biomass component, before (left bar, each mutant) and after (right bar, each mutant) the 
removal of latent pathways predicted by MOMA. (6) Percentage of elementary modes of each type disabled by latent pathway removal. Latent pathway removal 
renders a large number of elementary modes unavailable for use. Strikingly the majority of the eliminated modes involve biomass production. 
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Fig. S4: Range of possible growth rates when all optimally-inactive reactions are removed. The gray bars correspond to the range of realizable growth rates In each 
of the 52 knockout mutants when all reactions not active in the corresponding growth-maximizing state of the mutant are disabled. Purple circles correspond to the 
MOMA-predlcted growth rate In response to each knockout in this case. Orange circles correspond to the MOMA-predicted growth when only the latent pathways are 
removed, as presented in the main text. Although this extreme pathway removal is also predicted to Increase growth following a knockout perturbation, low-growth 
states remain available for the large majority of mutants. All growth rates are normalized by that of the optimal wild type. 
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Table S1 : Comparison between linear ROOM and integer ROOM. 



Integer ROOM 





+ 






Linear ROOM 








+ 


12 










12 


27 


1 














The rows and columns indicate the number of single-gene knockout perturbations for which the early post-perturbation 
growth-rate is larger {+), smaller (— ), or negligibly different (=) when the latent pathways are disabled beforehand 
compared to the case in which they are available. As in the main text, we define a change to be negligible if its 
magnitude is less than 1% of the wild-type growth rate. The majority of cases are on the diagonal, indicating that 
the two methods predict the same growth impact of latent pathway removal. There are no mutants for which the two 
methods predict a significant change in opposite directions. 



Table S2: Summary of the predicted impact of latent pathways under nutrient-specific regulatory constraints. 





MOMA 


ROOM 


Random 


Latent reactions for individual perturbations: 








All knockout perturbations 


+10.5 (13.2)% 


+1.5 (3.8)% 


+61.0 (9.4)% 


Significant differences * 


+18.4 (12.8)% 


+5.5 (5.9)% 


+61.0 (9.4)% 


Number of reactions removed 


267 (82) 


96 (45) 


902 (5) 


Simultaneously nonessential latent reactions: 
All knockout perturbations 
Significant differences ^ 
Number of reactions removed 


+8.8 (11.2)% 
+16.7 (10.4)% 
235 (76) 


+1.5 (3.8)% 
+5.8 (5.9)% 
82 (38) 





* By more than 1 % of the wild-type growth rate: 57% (MOMA), 26% (ROOM), and 1 00% (random) of the perturbations. 
tBy more than 1% of the wild-type growth rate: 52% (MOMA) and 24% (ROOM) of the perturbations. 



Corresponds to Table 1 of the main text for 46 knockout mutants in a modified /AF1 260 E. coli model in which 1 52 reactions have been disabled to 
account for known regulatory constraints in aerobic glucose medium conditions. 
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